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Background: Variable-number tandem repeat (VNTR) typing is a promising method to discriminate the 
Mycobacterium tuberculosis isolates in molecular epidemiology. The purpose of this study is to determine the optimal 
VNTR combinations for discriminating isolated M. tuberculosis strains in Korea. 

Methods: A total of 317 clinical isolates collected throughout Korea were genotyped by using the IS6110 restriction 
fragment length polymorphism (RFLP), and then analysed for the number of VNTR copies from 32 VNTR loci. 
Results: The results of discriminatory power according to diverse combinations were as follows: 25 clusters in 83 strains 
were yielded from the internationally standardized 15 VNTR loci (Hunter-Gaston discriminatory index [HGDI], 0.9958), 
25 clusters in 65 strains by using IS6110 RFLP (HGDI, 0.9977), 14 clusters in 32 strains in 12 hyper-variable VNTR loci 
(HGDI, 0.9995), 6 clusters in 13 strains in 32 VNTR loci (HDGI, 0.9998), and 7 clusters in 14 strains of both the 12 hyper- 
variable VNTR and IS6110 RFLP (HDGI, 0.9999). 

Conclusion: The combination of 12 hyper-variable VNTR typing can be an effective tool for genotyping Korean M. 
tuberculosis isolates where the Beijing strains are predominant. 

Keywords: Mycobacterium tuberculosis; Molecular Epidemiology; Minisatellite Repeats 



Introduction 

Tuberculosis (TB) is a contagious disease that develops 
from infection with Mycobacterium tuberculosis bacilli in 
droplets projected by coughing of active TB patients. DNA 
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typing of M. tuberculosis is an efficient tool to show the trans- 
mission link of TB 1 . In particular, IS6110 restriction fragment 
length polymorphism (RFLP) has been used for discrimina- 
tion of Korean M. tuberculosis isolates because of high IS6110 
copies and diverse patterns 2 4 . However, IS6110 RFLP has 
some disadvantages such as a lengthy turnaround time, dif- 
ficulty in comparison between laboratory DNA typing data, 
and uncertain discrimination of isolates containing similar 
fragments sizes for inter- or intra-strains 5 . 

Among more than 40 variable-number tandem repeat 
(VNTR) loci scattered on the M. tuberculosis chromosome 6 , 
15 and 24 VNTR loci have been proposed as the international 
standard 7 . However, the discriminatory power is not enough 
in countries that have a high proportion of Beijing type M. 
tuberculosis 812 . K strains that have 10 IS6110 RFLP copies 
are most dominant strain and often found in TB outbreak in 
Korea. K strains occupy about 4-5% out of any group of M. 
tuberculosis isolates isolated in Korea 3 ' 13 . K strains causes dif- 
ficulty in discerning the direct transmission link by IS 6110 
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RFLP typing only. 

Therefore, we attempted to optimize a combination of 
VNTR loci and IS6110 RFLP for discrimination of Korean M. 
tuberculosis isolates. 



Materials and Methods 

1. Strains 

A total of 317 strains were randomly selected among 2,400 
strains isolated from smear and culture-positive primary pul- 
monary TB patients that registered at the Public Health Cen- 
ter (PHC) of Korea between 2006 and 2011. When the strains 
were clustered in IS6110 RFLP or VNTR typing, we collected 



Table 1. Primer sequences and labeled dyes of 32 VNTR loci for multiplex PCR 



Multiplex 


Locus 


Alias 


Primer sequence 


Unit (bp) 


Dye 


Mixl 


2163a 


QUB-lla 


CGTGATGTTGATCGGGATGT/ACCCTGGAGTCTGGCATC 


69 


PET 




2996 


MIRU-26 


TAGGTCTACCGTCGAAATCTGTGAC/CATAGGCGACCAGGCGAATAG 


51 


VIC 




3192 


MIRU-31 


ACTGATTGGCTTCATACGGCTTTA/GTGCCGACGTGGTCTTGAT 


53 


NED 


Mix 2 


3820 




TGCGCGGTGAATGAGACG/ACCTTCATCCTTGGCGAC 


57 


PET 




2163b 


QUB-llb 


CCGATGTAGCCCGTGAAGA/AGGGTCTGATTGGCTACTCA 


69 


FAM 




4156 




ACCGCAAGGCTGATGATCC/GTGCATCTCGTCGACTTCC 


59 


NED 


Mix 3 


3232 




CCCCAGCCTTACGACTGA/GTCGGGCTTGGTGAAGG 


56 


FAM 




3336 




ATCCCCGCGGTACCCATC/GCCAGCGGTGTCGACTATCC 


59 


PET 




4052 


QUB-26 


GTGCCGGCCAGGTCCTTCC/CACCGCGTGTTTGACCCGAAC 


111 


NED 


Mix 4 


1955 


Mtub21 


AGACGTCAGATCCCAGTT/ACCCGACAACAAGCCCA 


57 


VIC 




4120 




GTTCACCGGAGCCAACC/GAGGTGGTTTCGTGGTCG 


57 


PET 




0424 


Mtub04 


CTTGGCCGGCATCAAGCGCATTATT/GGCAGCAGAGCCCGGGATTCTTC 


51 


FAM 


Mix 5 


0580 


MIRU-4 


GCGCGAGAGCCCGAACTGC/GCGCAGCAGAAACGTCAGC 


77 


FAM 




2165 


ETR-A 


AAATCGGTCCCATCACCTTCTTAT/CGAAGCCTGGGGTGCCCGCGATTT 


75 


NED 




1644 


MIRU-16 


TCGGTGATCGGGTCCAGTCCAAGTA/CCCGTCGTGCAGCCCTGGTAC 


53 


VIC 


Mix 6 


0960 


MIRU-10 


GTTCTTGACCAACTGCAGTCGTCC/GCCACCTTGGTGATCAGCTACCT 


53 


FAM 




3690 


Mtub39 


CGGTGGAGGCGATGAACGTCTTC/TAGAGCGGCACGGGGGAAAGCTTAG 


58 


VIC 




2074 


Mtub24 


TGTGTCACCTGACGATTTCAAGG/TGGCCGGCAAATAATGGATGC 


56 


PET 


Mix 7 


2347 


Mtub29 


GCCAGCCGCCGTGCATAAACCT/AGCCACCCGGTGTGCCTTGTATGAC 


57 


FAM 




3007 


MIRU-27 


TCGAAAGCCTCTGCGTGCCAGTAA/GCGATGTGAGCGTGCCACTCAA 


53 


NED 




2461 


ETR-B 


GCGAACACCAGGACAGCATCATG/GGCATGCCGGTGATCGAGTGG 


57 


VIC 


Mix 8 


2531 


MIRU-23 


CTGTCGATGGCCGCAACAAAACG/AGCTCAACGGGTTCGCCCTTTTGTC 


53 


VIC 




4348 


MIRU-39 


CGCATCGACAAACTGGAGCCAAAC/CGGAAACGTCTACGCCCCACACAT 


53 


NED 




3155 


QUB-15 


GCCAGCCGTAACCCGACCAG/GGGCCGGAAATTCGCAGTGG 


54 


PET 


Mix 9 


2687 


MIRU-24 


CGACCAAGATGTGCAGGAATACAT/GGGCGAGTTGAGCTCACAGAA 


54 


VIC 




0802 


MIRU-40 


GGGTTGCTGGATGACAACGTGT/GGGTGATCTCGGCGAAATCAGATA 


54 


NED 




2372 




ACCTCCGTTCCGATAATC/CAGCTTTCAGCCTCCACA 


57 


PET 


Mix 10 


2401 


Mtub30 


CGTCGTCGCCGAGCTGGATT/CACCGGGGCTGGCAGCTAAG 


58 


PET 




3171 


Mtub34 


GGTGCGCACCTGCTCCAGATAA/GGCTCTCATTGCTGGAGGGTTGTAC 


54 


NED 




0577 


ETR-C 


GTGAGTCGCTGCAGAACCTGCAG/GGCGTCTTGACCTCCACGAGTG 


58 


VIC 


Single 


0154 


MIRU-2 


TGGACTTGCAGCAATGGACCAACT/TACTCGGACGCCGGCTCAAAAT 


53 


FAM 


Single 


2059 


MIRU-20 


TCGGAGAGATGCCCTTCGAGTTAG/GGAGACCGCGACCAGGTACTTGTA 


77 


FAM 



VNTR: variable-number tandem repeat; PCR: polymerase chain reaction. 



60 



Tuberc Respir Dis 2014;76:59-65 



www.e-trd.org 



VNTR typing of Korean M. tuberculosis isolates 



Tuberculosis and ^% 
Respiratory Diseases I w\ 



the epidemiological information of the strains from the PHCs 
and each person who had the clustered strain by documents 
and telephone calls. 

2. IS6II0DNA fingerprinting 

For all 317 isolates, DNA isolation and IS6110 RFLP typing 
were performed as described previously 2 . An RFLP cluster 
was defined by completely identical patterns after analysis 
of BioNumerics version 5.1 software (Applied Maths, Kortrjk, 
Belgium) within two or more isolates. K strains were identified 
according to the previously reported definition 4 13 . 

3. Variable number of tandem repeats 

We selected 32 VNTR loci for this study 7 9 . Primer sets for 
polymerase chain reaction (PCR) were prepared for the 32 
VNTR loci. PCR was conducted as described previously 14 16 
and VNTR 4052, 3155, and 2074 were designated by the 
Research Institute of Tuberculosis in Japan (Table 1). Each 
primer set was labelled with 4 kinds of fluorescent dye for 
capillary sequencer analysis. The capillary sequencer, a 3500 
Genetic analyzer (Applied Biosystems, Foster City, CA, USA), 
was used for measuring the precise size of PCR products for 
each VNTR loci. However, in cases of large fragment sizes over 
12,000 bp or an ambiguous size in the capillary sequencer, we 
measured the size by 1% agarose gel electrophoresis. Clusters 
were found as the result of comparison of genotyping data by 
BioNumerics version 5.1 software (Applied Maths). For PCR 
amplification, we mainly used Ex-Taq polymerase (Takara, 
Tokyo, Japan), except for Mix4, Mix7, and Mix9, which used 
KOD FX Taq polymerase (TOYOBO, Tokyo, Japan). PCR 
conditions were as follows: pre-denaturation for 5 minutes 
at 94°C, denaturation for 30 seconds at 94°C, annealing for 
30 seconds at 63°C, extension for 1 minute at 72°C, and post- 
extension for 10 minutes at 72°C, except for Mix9 that had an 
annealing temperature of 60°C. 

4. Comparison of discriminatory power 

The allelic diversity {h) of the 32 VNTR loci was calculated 
by the following formula: h=l -5X 2 , where x 4 is the frequency of 
the ith allele at the locus 17 . The discriminatory power of each 
method was calculated by the Hunter-Gaston discriminatory 
index (HGDI) 18 . 

5. Ethical considerations 

Ethical clearance was obtained from the Ethics Review 
Committee of Korean Institute of Tuberculosis. 



Results 

1 . Analysis of the 32 VNTR loci 

As a result of VNTR typing of the 317 M. tuberculosis strains 
isolated from Korea, VNTR loci showing a high h value of over 
0.6 were VNTR 3232, 3820, 4120, 3336, 2163b, 0424, 1955, 
4052, 3192, 4156, 2996, and 2163a. Interestingly, VNTR 0802, 
3690, 2165, and 0960 revealed a high h value (over 0.6) only in 
non-Beijing M. tuberculosis isolates (Table 2). 

Two PCR products amplified from VNTR 3802 and 0580 
(MIRU04) were an unusual size, and were regarded as un- 
countable repeated numbers (Table 2). The imprecise PCR 
product of the VNTR 3802 loci was located between 8 and 9 
copies, and that of VNTR 0580 was located between 3 and 4 
copies of the repeated unit. Four PCR products in VNTR 3820, 
4052 (QUB26), and 2372 had multiple bands. Non-amplified 
PCR products were found in VNTR 3232, 4120, and 2163b 
(QUBllb), 3192 (MIRU31), 2163a (QUBlla), 2165 (ETR-A), 
and 2074 (Mtub24). 

2. Epidemiological linkage of the VNTR clusters 

Six clusters in 13 patients were found after VNTR typing of 
the 32 loci, and only 3 clusters in 6 patients were found when 
additionally analysed by IS6110 RFLP typing (Table 3). Out 
of the 3 clusters, we found only one definite epidemiological 
linkage through personal contact. Patients with the V32C6 
cluster were brother and sister. Patient 06-1731 with the clus- 
ter was the brother who developed TB in 2006, and was cured 
completely in 2007. Patient 08-148, who was the sister of pa- 
tient 06-1731, had developed TB, and was cured in 2008. Even 
though random selection of PHC strains, we found the strains 
isolated from brother and sister TB-developed in different 
year by accident. 

Two patients with V32C3 and IS6110 clusters had lived in 
Jeju province, even though they were not acquainted with 
each other. The two strains (06-890 and 09-1184) with the 
C32C4 cluster and identical IS6110 RFLP types were the K 
strain that is the most frequent endemic strain in Korea (4-5% 
in any TB population). The other strain (10-335) with the 
C32C4 cluster was the K family that exhibited a difference of 
only one band. Interestingly, another patient (10-179) with the 
V32C1 cluster also had the K strain. However, we found dif- 
ferent copies in VNTR 2163a, 3232, and 3820 compared with 
that in V32C4 (data not shown). 

3. Discriminatory power of VNTR and 1S6110 RFLP typ- 
ing 

When the internationally standardized 15 VNTR loci were 
applied to the 317 strains, we found 25 clusters in 83 strains 
(HGDI, 0.9958) (Table 4), and at least 7 VNTR loci (3232, 
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Table 3. Characteristics of 6 clusters isolated from VNTR typing analysis in 32 VNTR loci 



VNTR 
cluster 


Year 


Region 


ID 


No. of IS61 10 
copies 


IS61 1 0 RFLP pattern 




V32C1 


2006 


Ansan 


06-681 




9 










2010 


Jungsun 


10-179 




10 




mil ii ii 




V32C2 


2007 


Pyongtaek 


07-230 




12 




mi m ■■ 






2007 


Pyongtaek 


07-977 




13 




nil ni i n 




V32C3 


2006 


Jeju 


06-481 




17 




ikiiii ii i tm 






2006 


Jeju 


06-571 




17 




Ml I ILI II 1 mM 




V32C4 


2006 


Ansan 


06-890 




10 










2009 


Gosung 


09-1184 




10 




MINI II II 

II 1 Ml II II 






2010 


Pohang 


10-335 




11 




HHH II II 




V32C5 


2007 


Anyang 


07-2256 




18 




mil ii dm inn i 






2007 


Tongyoung 


07-1645 




19 




mi 1 1 1 d in iiiii i 




V32C6 


2006 


Busan 


06-1731 




16 




i ii i u 1 1 mi: 






2008 


Busan 


08-148 




16 




1 11 1 1 1 1 1 1 I Ml 




VNTR: variable-number tandem repeat; RFLP: restriction fragment length polymorphism. 






Table 4. Comparison of discriminatory power according to DNA typing method alone or combination 








No. of 
DNA types 


No. of 
unique types 




No. of 
clusters 


Total isolates Maximum isolates 




Typing method 




belonging to belonging to a 
clusters cluster 


HGDI 


IS6ii0RFLP 




277 




252 




25 


65 13 


0.9977 


12 VNTR 




299 




285 




14 


32 4 


0.9995 


15 VNTR loci 




257 




232 




25 


83 16 


0.9958 


32 VNTR 




310 




304 




6 


13 3 


0.9998 


32VNTR+IS6iiORFLP 


314 




311 




3 


6 2 


0.9999 


12VNTR+IS6ii0RFLP 


310 




303 




7 


14 2 


0.9999 



Fifteen variable-number tandem repeat (VNTR) loci are of international standards, and 12 VNTR are hyper-variables. 
RFLP: restriction fragment length polymorphism. 



3336, 3820, 4120, 2372, 4348, and 3155) were needed to be 
equivalent to the discriminatory power of the 32 VNTR loci 
(data not shown). However, in cases of the hyper-variable 12 
VNTR loci (3232, 3820, 4120, 3336, 2163b, 0424, 1955, 4052, 
3192, 4156, 2996, and 2163a), 14 clusters in 32 strains (HGDI, 
0.9995) were found, and we needed only 4 VNTR loci (4348, 
2165, 0577, and 2372) to be equal to the discriminatory power 
of the 32 VNTR loci (data not shown). Furthermore, in the 
IS6110 RFLP typing, 25 clusters in 65 strains (HGDI, 0.9977) 
(Table 4) were obtained, and we needed at least 14 VNTR loci 
(4120, 3820, 2372, 1955, 0424, 4052, 3155, 2163b, 2163a, 0577, 
0802, 4348, 3192, and 2996) for 3 clusters as a result of both 
the 32 VNTR loci and IS6110 RFLP typing (data not shown). 
When the 12 hyper-variable VNTR loci were applied to the 
317 strains, followed by IS6110 RFLP typing, the discrimina- 



tory power was almost identical to the result of both the 32 
VNTRloci and IS6ii0RFLP typing (Table 4). 



Discussion 

VNTR typing has been proposed as an alternative to IS6110 
RFLP typing that bears intrinsic drawbacks for M. tuberculo- 
sis DNA typing. However, the disadvantage of VNTR typing 
is low discriminatory power with few combinations of VNTR 
loci. Therefore, we tried to find an optimal combination of 
DNA typing methods to discriminate M. tuberculosis isolated 
from Korea, which is accompanied by a high proportion of 
Beijing type 3 4 . Furthermore, compared with the surrounding 
countries, there is a higher proportion of RD181 among the 
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ancient Beijing types 4 . These peculiar characteristics of Ko- 
rean M. tuberculosis strains lead to an unclear distinction of M. 
tuberculosis isolates with only the 15 international standard 
VNTRloci. 

Recently, the utility of hyper-variable VNTR loci has been 
used to compensate for the lack of discriminatory power of 
the 15 VNTR loci 10 ' 11 ' 16 . Hyper-variable VNTR loci, including 
3232, 3820, 4120, 3336, and 2163a, also revealed a high h value 
in this study, which was similar to that in other studies 9 " 11,16 . In 
the report Murase et al. 9 , more than 4% of theses 5 loci had 15 
or more copies, leading to difficulty in interpreting the exact 
copy number. We also found that 2-4% had 15 or more copies 
in VNTR 2163b, 3336, and 4120, and 13-16% had 15 or more 
copies in VNTR 3820, and 3232. These hyper-variable VNTR 
loci are excluded in the international standard 15 VNTR loci 7 
and JATA 12 9 because of the absence of PCR products, PCR 
products that are difficult to interpret with 15 or more cop- 
ies, and amplification of multiple alleles. Non-amplification 
of PCR products not only occurred for hyper-variable VNTR 
loci but also general VNTR loci such as 2165 (ETR-A), and 
2074 (Mtub24) in this study (Table 2). The major problem of 
hyper-variable VNTR loci was a high copy number of the re- 
peats, which required additional analysis such as agarose gel 
electrophoresis. However, the copies in these hyper-variable 
VNTR loci are so diverse that they have a higher discrimina- 
tory power that is too valuable to exclude. Iwamoto et al. 16 
also recommended that these hyper-variable VNTR loci for 
second-line typing of clustering following the international 
standard 15 loci. 

MIRU 40, Mtub 21, VNTR 4156, Mtub 04, QUB26, and 
QUBllb also revealed high h values in a previous study of 
South Korean M. tuberculosis 19 . In this study, we obtained a 
good result using hyper-variable VNTR loci, and IS6110 RFLP 
was useful as a secondary tool to discriminate the clusters. 

An intriguing characteristic was found after comparing the 
mode of VNTR loci among Korea, Taiwan, and Japan 9 20 . Be- 
tween Taiwan and Korea, there were some differences in the 
mode of copies of VNTR loci 0154, 3192, 2163b, 4052, 0424, 
1955, 2347, and 2401. In particular, 5 copies of VNTR 0154, 
0424, and 1955, and 3 copies of VNTR 2347 were modes in 
Taiwan but rare in Korea. Between Japan and Korea, there 
were some differences in VNTR 2163b, 1955, 3155, 3232, 
3820, and 4120. Notably, 5 copies of VNTR 3155 and 4 copies 
of VNTR 3336 were more frequent in Japan but rare in Korea. 
These differences may be a clue for differentiation among the 
three countries. Compared with Beijing strains, we found that 
VNTR 2165, 0960, and 2074 loci had excellent higher h values 
for the non-Beijing strains, indicating that these VNTR loci 
may be useful to discriminate M. tuberculosis strains in coun- 
tries with a high proportion of non-Beijing strains. The optimal 
combination of VNTR loci may be different depending on the 
proportion of Beijing strains or non-Beijing strains in each 
country. 



We could not find any cluster consisted of only multi-drug 
resistant (MDR) strains in this study. The strains included in 
this study were selected randomly among strains collected 
from PHCs. Therefore, most of them (281 strains) were pan- 
susceptible, only a few of strains (13 strains) were MDR, and 
the remains (23 strains) were any drug resistant with non- 
MDR. 

Even though we analysed 32 VNTR loci for discrimination 
of 317 Korean M. tuberculosis isolates, we found 6 clusters. 
Among them, 3 clusters were not found to be clusters when 
additionally analysed by the IS6110 RFLP typing method, 
indicating that IS6110 RFLP typing is very useful for sub-clas- 
sifying VNTR clusters. In terms of cost efficiency, it is difficult 
to use 32 VNTR loci for discrimination of M. tuberculosis, and 
the 15 international standard VNTR loci do not have satisfac- 
tory discrimination power for Korean strains. Inevitably, we 
need hyper-variable VNTR loci and the additional IS6110 
RFLP typing method for effective discrimination of Korean M. 
tuberculosis strains. 

The combination of 12 hyper-variable VNTR typing can be 
an effective tool for genotyping Korean M. tuberculosis iso- 
lates in which Beijing strains are predominant. 
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